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ABSTRACT 



43 . 

Aims. This paper reports on H-band interferometric observations of Betelgeuse made at the three-telescope interferometer IOTA. We 
image Betelgeuse and its asymmetries to understand the spatial variation of the photosphere, including its diameter, limb darkening, 
O 1 effective temperature, surrounding brightness, and bright (or dark) star spots. 

Methods. We used different theoretical simulations of the photosphere and dusty environment to model the visibility data. We made 
C/} , images with parametric modeling and two image reconstruction algorithms: MIRA and WISARD. 

Results. We measure an average limb-darkened diameter of 44.28 + 0.15 mas with linear and quadratic models and a Rosseland 
diameter of 45.03 + 0.12 mas with a MARCS model. These measurements lead us to derive an updated effective temperature of 
3600 + 66 K. We detect a fully-resolved environment to which the silicate dust shell is likely to contribute. By using two imaging 
reconstruction algorithms, we unveiled two bright spots on the surface of Betelgeuse. One spot has a diameter of about 11 mas and 
accounts for about 8.5% of the total flux. The second one is unresolved (diameter < 9 mas) with 4.5% of the total flux. 
Conclusions. Resolved images of Betelgeuse in the H band are asymmetric at the level of a few percent. The MOLsphere is not 
detected in this wavelength range. The amount of measured limb-darkening is in good agreement with model predictions. The two 
spots imaged at the surface of the star are potential signatures of convective cells. 
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■ 1. Introduction the size for SI and S2. An outer radius of -4.5 and ~7 arc- 
sec were derived, respectively, although the latter is inconsistent 

Located in the Orion constellation, Betelgeuse (a Orionis) is a with other meaS urements made by the CARMA interferometer. 
K> , red supergiant (hereafter RSG) of spectral type M2Iab. It is one 
^ ■ of the brightest stars at optical wavelengths and has the second 

C3 | biggest angular diameter (~43 mas, Perrin et al., 2004) after 1.2. A dusty environment 
R Doradus (Bedding et al., 1997). Classified as semi-regular, 

it shows periodicity in its brightness changes, accompanied or A shell of dust was first detected with heterodyne interferom- 

sometimes interrupted by various irregularities (Guinan, 1984). etry at 11 //m by Sutton et al. (1977). Quite a few estimates 

According to a recent reanalysis of a Hipparcos satellite dataset, of the inner radius of Betelgeuse's dust shell have since been 

its distance is now estimated at 197 ±45 pc (Harper etal, 2008). made. Bester et al. (1991) used a model with an inner ra- 

The observations up to now have identified at least 7 components dius of 0.9arcsec (-45 R*) to explain both their llyum het- 

of the complex Betelgeuse atmosphere: two outer shells, a dust erodyne interferometry (ISI) and older speckle observations by 

environment, a chromosphere, a gaseous envelope, a molecular Sutton et al. (1977) and Howell et al. (1981). In their spatially- 

shell also known as MOLsphere (Tsuji, 2000a), and finally the resolved mid-infrared (mid-IR) slit spectroscopy, Sloan et al. 

photosphere. Some of these components are not symmetric, and (1993) find no silicate emission within the central arcsecond 

some are overlapping. around Betelgeuse 1 . Danchi et al. (1994) find from 11.15/zm 

ISI data that the inner radius must be 1.00 ± 0.05arcsec, i.e. 

roughly 50 R*. But this result disagrees with later findings by 

1.1. Two outer shells Skinner et al. ( 1 997), who claim an inner radius of not more than 

rr u ti /ci a oo\ u u u a • u 0.5 arcsec. This environment of dust is also reported by Tatebe 

Two shells (SI and S2) have been observed m absorption in v j 

12 C 16 and 13 C 16 at 4.6 /mi by Bernat et al. (1979) but their 

spatial extent was not directly determined. Phoenix 4.6 yum spec- 1 They also remark that this actually argues against a spherical distri- 

tra obtained by Harper et al. (2009) led to an estimation of bution of the dust. 
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et al. (2007) as the origin of a visibility drop at short baselines 
of more than 40% in new ISI measurements. 

1.3. Chromosphere 

Lambert & Snell (1975) interpreted Betelgeuse spectra and 
found that the infrared excess could be explained by a chro- 
mospheric emission. More recently, the ultraviolet (UV) ob- 
servations of the Hubble Space Telescope revealed a spatially 
resolved chromosphere measuring twice the size of the opti- 
cal photosphere in the continuum (Gilliland & Dupree, 1996; 
Uitenbroeket al., 1998). 

1.4. The asymmetric gaseous envelope 

The outer envelope of the O-rich Betelgeuse (0/C=2.5, Lambert 
et al., 1984) has been the subject of extensive studies. VLA ob- 
servations revealed an extended, cool (-1000-3000 K) and irreg- 
ularly shaped gaseous atmosphere of several stellar radii which 
coexists with the chromospheric gas (Lim et al., 1998). More re- 
cently Kervella et al. (2009) found evidence for the presence of 
a bright plume, possibly containing CN, extending in the south- 
west direction up to a minimum distance of six photospheric 
radii. Such a plume could be formed either by an enhanced mass 
loss above a large upwards moving convective cell or by rota- 
tion, through the presence of a hot spot at the location of the 
polar cap of the star. 

1.5. The MOLsphere 

With the Stratoscope II experiment, strong water-vapor bands 
have been detected in the spectrum of Betelgeuse (Danielson et 
al., 1965) but at that time no explanation on their origin was 
given. More recently, Tsuji (Tsuji, 2000a,b) proposed the exis- 
tence of a molecular layer at 0.45 above the photosphere, a 
MOLsphere, of effective temperature 1500 ± 500 K and com- 
posed of water vapor and CO. Using infrared interferometry, 
Perrin et al. (2004) showed that K, L and 11.15 yum data could 
be simultaneously modeled with a gaseous shell of temperature 
2055 ± 25 K located 0.33 7?* above the photosphere while Ryde 
et al. (2006) locate the water vapor inside the photosphere by an- 
alyzing high-resolution spectra. A model based on mid-infrared 
measurements showed that the MOLsphere could contain not 
only water vapor but also SiO and AI2O3 (Verhoelst et al., 2006; 
Perrin et al., 2007). This latter component is a type of dust con- 
densing at high temperature and on which SiO can be adsorbed. 
It is a nucleation site at high temperature for silicate dust. This 
is a possible scenario to explain dust formation and to connect 
molecular layers and dusty environments. 

1.6. The convective photosphere 

Perrin et al. (2004) modeled the photosphere of Betelgeuse as a 
uniform disk of 3641 ± 53 K temperature with a 42 mas diame- 
ter. For some wavelengths, the diameter can vary with time. At 
11.15yum, Townes et al. (2009) measure a systematic decrease 
of the diameter of Betelgeuse during the period 1993-2009 that 
might point toward long-term dynamics at its surface. Moreover, 
the more recent hydro-radiative simulations of RSGs (Freytag et 
al., 2002) foresee the presence of convection cells whose size is 
comparable to the stellar radius, as suggested by Schwarzschild 
(1975). 



1.7. The spots observed on the surface of Betelgeuse 

Thanks to the pupil masking technique, Wilson et al. (1992, 
1997) and Buscher et al. (1990) have detected the presence of 
spots on the surface of Betelgeuse that represent up to 15-20% 
of the observed flux in the visible. The most likely origin of 
these asymmetries is a signature of the convection phenomenon. 
However, the hypotheses of inhomogeneities in the molecular 
layers or the presence of a transiting companion are not ruled 
out (Buscher et al., 1990). These observed spots are few, 2 or 
3, and their characteristrics change with time (Wilson et al., 
1992). Confirming previous results, Young et al. (2000) observed 
spots in the visible but show a fully centra- symmetric picture of 
Betelgeuse in the near infrared (at 1.29 yum). Lately Tatebe et 
al. (2007) also report the observation of an asymmetry at 11.15 
jim located on the southern edge of the disk of Betelgeuse. In 
order to deepen our knowledge on RSGs, it is critical to observe 
spots at high angular resolution. Many unknowns remain such as 
their size, location, chemical composition, dynamical properties, 
lifetime and origin. 

To investigate these questions, we report in this paper on in- 
terferometric observations obtained with the IONIC 3 -telescope 
beam combiner at the IOTA interferometer in October 2005. 
Sect. 2 describes the observations and the data reduction. In 
Sect. 3, the fully resolved environment is discussed. The physi- 
cal diameter and limb darkening of Betelgeuse are investigated 
in Sect. 4 . In Sect. 5, images of Betelgeuse reconstructed with 
the MIRA and WIS ARD algorithms are compared to a paramet- 
ric fit. Conclusions are gathered in Sect. 6. 



2. Observation and Data reduction 

The data were acquired with the 3 -telescope interferometer 
IOTA (Infrared Optical Telescope Array, Traub et al., 2003) lo- 
cated on Mount Hopkins in Arizona. This instrument has an 
L shape whose arms are oriented northeast and southeast with 
lengths of 35 and 15 meters respectively. In order to sample the 
uv plane, the three telescopes are movable on these arms and 
can be locked at 17 different locations. The maximum projected 
baseline is thus 38 meters corresponding to a resolution 2 of 9 
mas in H band. Light coming from three apertures (siderostats of 
0.45 m diameter) is spatially filtered by single-mode fibers which 
clean the wavefront and remove atmospheric corrugations that 
affect the contrast (Coude du Foresto et al., 1998; Perrin et al., 
1998). The fringes are scanned with piezo mirrors. The beams 
are then combined with IONIC (Berger et al., 2003). This inte- 
grated optics (IO) component combines 3 input beams in a pair- 
wise manner. The combination is coaxial. Light was split in two 
orthogonal polarization axes with a Wollaston prism. Finally the 
fringes were imaged with an acquistion camera using a PICNIC 
detector (Pedretti et al., 2004). We observed a Ori in the H 
narrow-band filter (Ao=l.64 yum, AA =0.10 yum) during 6 days 
using 5 different configurations of the interferometer in order 
to get a rich uv coverage (see Table 1 and Fig. 1). To calibrate 
the instrumental transfer function, we observed alternatively a 
Orionis and a calibrator (HD 36167). 

Data reduction was carried out using an IDL pipeline 
(Monnier et al., 2004; Zhao et al, 2007). For the visibilities, 



2 Here, we define the resolution as | which is given by the cutoff 
spatial frequency of the Optical Transfer Function (OTF), A being the 
observing wavelength and B the projected baseline. We mention that 
alternative values exist such as the FWHM of the central fringe. 
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Date 


Configuration 


Maximum projected baseline 


Resolution associated with the 






in meters 


maximum projected baseline in mas 


WLIUUCI /, ZUUJ 


A< R< cc\ 


U.J 


S1 7 
J 1 . / 


October 8, 2005 


A5-B15-C0 


15.2 


22.1 


October 10, 2005 


A20-B15-C0 


19.5 


17.2 


October 11,2005 


A20-B15-C0 


19.6 


17.1 


October 12, 2005 


A25-B15-C0 


27.2 


12.4 


October 16, 2005 


A30-B15-C15 


31.2 


10.8 



Table 1. Log of the observations. Configurations represent the station locations of the telescopes A,B and C along the northeastern arm (A and C) 
and southeastern arm (B) of IOTA. 
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Fig. 1. uv coverage of the Betelgeuse observations. 



we measured the power spectrum of each interferogram (propor- 
tional to the target squared visibility) after correcting for inten- 
sity fluctuations and subtracting out bias terms from read noise, 
residual intensity fluctuations, and photon noise (Perrin, 2003). 
Then, the data pipeline applies a correction for the variable flux 
ratios for each baseline by using a flux transfer matrix (Coude du 
Foresto et aL, 1997; Monnier, 2001). Finally, raw squared visi- 
bilities are calibrated using the raw visibilities obtained by the 
same means on the calibrator stars. 

Closure phases are computed as the phase of the complex 
triple amplitude, also called bispectrum (Baldwin et aL, 1996). 
These measurements were recorded when interferences occurred 
on the three baselines taking advantage of the baseline boot- 
strapping technique (Pedretti et aL, 2005). This method makes 
it possible to record fringes on a long baseline at low visibil- 
ity when fringes on the shorter baselines are detected. Pairwise 
combiners (such as IONIC3) can have a large instrumental offset 
for the closure phase. This is usually calibrated by using a point- 
source calibrator. The instrumental closure phase of IONIC3 
drifted by less than 1 degree over many hours, thanks to the com- 
pact size of the IO component. For both the squared visibilities 
and the closure phase, the errors were calculated with the boot- 
strap resampling method. This technique involves repeated re- 
estimation of a statistical parameter using random samples with 
replacement from the original data. 

IO combiners always have some low level of photon 
crosstalk. Here, the crosstalk phenomenon is defined as the un- 
desired transmission of flux (coherent or not) coming from one 
telescope to a combiner (therefore baseline) involving two other 
telescopes. This is a potential source of bias, especially when 



measuring small visibilities, that needs to be carefully calibrated. 
We investigated whether the level of visibilities at high frequen- 
cies is due to a crosstalk contamination from a shorter baseline. 
From past crosstalk measurements, we simulated a pessimistic 
crosstalk where the high contrast obtained on a short baseline 
contaminates a lower contrast measured on a longer baseline. 
The results of these simulations show that the crosstalk cannot 
account for more than half a percent of the high spatial frequency 
visibilities we measured, i.e. the maximum potential visibility 
bias is within error bars. Other observations of Arcturus with 
IONIC/IOTA (Lacour et aL, 2008) revealed a minimum level 
of 1 % of visibility at high spatial frequency. With a contrast up 
to -4.5%, the visibilities at high spatial frequency observed on 
Betelgeuse cannot be explained by a crosstalk effect and we fi- 
nally conclude they undoubtedly are of astrophysical origin. 

3. Detection of the circumstellar environment of 
Betelgeuse 

The whole data set comprizes 393 squared visibility and 131 
closure phase calibrated measurements. We tried several ap- 
proaches in order to fit the data. We first compare the data with 
a simple disk model and then discuss the circumstellar environ- 
ment and limb darkening. 

3.1. Uniform disk and a fully resolved environment 

A first model that can be fit to interferometric data is the uniform 
disk (UD). It has only one parameter, the UD diameter cp UD . It 
is used to give a first approximation to the object's angular size. 
Only first lobe data are used. Higher order lobes contain infor- 
mation on limb darkening and on smaller details. A loss of vis- 
ibiity at low spatial frequencies (i.e. at short baselines) needs to 
be modeled by scaling the visibility function. Thus the squared 
visibility is 

,v (B ,«p=Lxr ¥ (H!swf (1) 

I J Ao -f \ n<f>BJA ) I 

where B ± = Vw 2 + v 2 and J\ is the second Bessel function of the 
first kind. The constant factor co represents the loss of contrast at 
null baseline due to the presence of a fully resolved environ- 
ment, i.e., circumstellar light, from a region outside the stellar 
disk. Consequently, the visibility of the central object at the null 
baseline is not equal to one anymore but to this scaling factor. 
The fit is shown in Fig. 2. 

The second part of the first lobe dataset is well reproduced 
by this model. Consequently, the reduced x 2 changes from 42 
(for one UD) to 16 (UD + environment). This model gives a UD 
diameter of 42.41 ± 0.07 mas and a coefficient co of 0.96 ± 0.01, 
meaning that the environment is responsible for about 4 % of the 
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3.2. Extended emission 

Extended broadband emission at a scale of -100 mas and more 
has been detected around Betelgeuse in mid-IR, radio and at UV 
wavelengths. The mid-IR is caused by thermal dust emission 
(see e.g. Danchi et al., 1994). The radio emission mostly arises 
from free-free processes in the gaseous envelope (Harper et al., 
2001). The UV emission is a result of a combination of scattered 
and fluorescent line photons, electron collisional excitation, and 
bound-free recombination in the chromosphere and/or ionized 
wind (Harper et al, 2001; Lobel & Dupree, 2001). Although the 
H band covers an opacity minimum of the combined free-free 
and bound-free emission from both chromosphere and (partly) 
ionized wind in the gaseous envelope, we can not exclude these 
emission processes as the origin of the observed environment 
without detailed calculations based on the models which were 
constructed to explain the UV and radio observations (Harper 
et al., 2001, e.g.). However, as is shown in this section, scattered 
light on the inner edge of the dust shell provides a natural expla- 
nation of the observed loss of fringe contrast. The calculation of 
a possible gaseous emission component at these wavelengths is 
left to a future paper. 
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Fig. 2. Top: Fitting of a UD plus an environment to the first lobe of 
Betelgeuse. Bottom: Enlargements of the two parts of the first lobe. 
Squared visibility at null baseline is equal to 0.93 (meaning 0.96 in 
visibility). 



total flux. That UD diameter is about one mas smaller than the 
previous measurement of 43.26 mas made by Perrin et al. (2004) 
based on data in K band. 

A loss of constrast can also exist when one observes an ob- 
ject so bright that the detector works in a non-linear regime or 
even saturates. This case was possible given the high flux of 
Betelgeuse in H band (m#=-3). Although we used diaphragms 
to decrease the stellar flux when we acquired the data, it is con- 
ceivable that we approached the non-linear regime of the detec- 
tor due to variability in the observing conditions. After a visual 
analysis of the interferograms and the study of a photometric 
estimator, we concluded that our data do not suffer from non- 
linearity or saturation from the detector. Moreover, for all the 
observed stars during that campaign, no other has shown a simi- 
lar effect. Thus, the effect does not have an instrumental origin. 

The necessity to account for a visibility scaling factor is 
hence the first evidence of the presence of a fully resolved en- 
vironment in the near infrared. Even the shortest baseline data 
do not appear to be contaminated by correlated flux from this 
extended component at the level of the visibility accuracy. A 
rough estimate of the minimum size of this environment can be 
established by considering no departure from the uniform disk 
model is detected at short baseline. Assuming simple models for 
the environment like a Gaussian distribution or a uniform disk 
distribution and a maximum contribution of 1 % to the visibility 
residuals, minimum sizes of around 100 and 400 mas are respec- 
tively derived. The size of the environment is therefore probably 
at least a few hundred mas. The visibility drop could stem from 
an extended environment in the upper atmosphere or to a dusty 
shell. We study these hypotheses in the following. 



3.3. Dust shell model and scattered light contribution 

The dust shell was modeled using the proprietary spherical ra- 
diative transfer code modust (Bouwman et al., 2000; Bouwman, 
2001). Under the constraint of radiative equilibrium, this code 
solves the monochromatic radiative transfer equation from 
UV/optical to millimeter wavelengths using a Feautrier type so- 
lution method (Feautrier, 1964; Mihalas, 1978). The code per- 
mits several different dust components of various grain sizes 
and shapes, each with their own temperature distribution. In 
Betelgeuse, the main spectral dust feature is that of amorphous 
Magnesium-Iron silicates. Input parameters were derived by 
modeling the observed IR spectrum in Verhoelst et al. (2006). 
This kind of modeling of optically-thin dust shells, which is 
based on (a limited part of) the IR spectrum and on the opti- 
cal obscuration, is slightly degenerate between inner radius and 
mass-loss rate, especially when the exact dust composition is not 
known. 

Regardless of the exact inner radius, we find that any dust shell 
model explaining the observed silicate feature at mir-IR wave- 
lengths produces roughly 4% scattered stellar light in the H band. 
The crucial question is: which fraction of this scattered light 
makes it into the IONIC fibers (with a field of view of 750 mas 
diameter) ? This depends solely on the assumed inner radius. 
Fig. 3 shows the amount of excess making it into the fiber for 
a model with a 13R* inner radius. Squared and Gaussian pro- 
files have been choosen to simulate the fiber field of view. We 
used diaphragms to avoid saturation of the detector during the 
observations. The diaphragms were located just before the fiber 
injection units. We thus take into account that the width of the 
fiber field of view was changing depending on the size of the 
diaphragm. For the short baselines where the fully resolved flux 
has been detected, the same diaphragm was employed. So we 
simulate our profiles for two realistic widths of the fiber field 
of view, with and without diaphragm: respectively 2.85 and 0.76 
arcseconds. In both cases, we expect a roughly 2% fully-resolved 
flux. Our observations may be compatible with a smaller inner 
radius, resulting in more scattered light making it into the fiber. 
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Fig. 3. Amount of beam integrated flux and dust shell intensity com- 
pared to the fiber profile. Black line: beam integrated flux, normalized 
to the photospheric flux, as a function of beam radius, either in units of 
dust shell inner radius (bottom axis) or in arcseconds assuming a 13 R* 
dust shell inner radius (top axis). Grey line: the intensity profile of the 
dust shell, normalized to the central photospheric intensity. The dashed 
line represents the Gaussian fiber profile. 

4. Limb-darkening and physical diameter 

The limb darkening effect is a gradual decrease of surface bright- 
ness from the center to the stellar limb caused by the decrease 
of temperature with altitude: the optical depth of 1 is reached 
for lower altitudes at the center of the stellar disk than at the 
edge. Taking this effect into account is mandatory for an accu- 
rate measurement of the diameter of the photosphere. To do so, 
we used parametric models (see Claret, 2004; Hestroffer, 1997) 
and a physical model derived from a MARCS simulation of the 
Betelgeuse atmosphere. 

4.1. Parametric models 

In order to compare different observations, one needs standard 
models of the stellar intensity distribution. The simplest ones 
are the linear and quadratic limb darkening laws. The visibility 
of a point- symmetric object is a Hankel transform (whose ker- 
nel is a Bessel function) normalized by the total spatial intensity 
distribution. For the monochromatic case, the visibility is 



poo B j_ 

J I(r,A)J (27T—r)rdr 



(2) 



tf °°I(r,A)rdr\ 



where B ± is the projected baseline of the interferometer, r the 
distance to the star center, I(r, A) the spatial intensity distribution 
at the observing wavelength A and Jq is the first Bessel function 
of the first kind. We present now the calculation of the visibil- 
ity for a quadratically limb darkened disk. For an achromatic 
quadratic limb darkening, the spatial intensity distribution is a 
function of jj, the cosine of the angle between the radius vector 
and the line of sight. It can be written as follows: 



m=l-A(l-fi)-B(l-fi) 2 



(3) 



with \i = Si - where (pLD is the limb darkened disk 

diameter. With B=0, one gets the linear limb darkening case. 



Substituting the intensity equation in the visibility expression, 
the Mellin transform appears. For the general case of a power 
law of index S, = \a 5 , this can be written as follows 

(Hestroffer, 1997): 



W(jc) = 2v 



I ft 1 -'' 
Jo 



2 ) v 1 Jo(x/j)fid/2 



= r(v + i)x 



\Jv(x)\ 

(x/iy 



(4) 



where J v is the Bessel function of the first kind of index v, 
defined as v = 6/2 + 1. 

The expression for the visibility of a quadratic limb- 
darkening is thus finally: 



|/ife)l 



^ \J3/2(Z)\ \J 2 (Z)\ 



VAz) 



with a 



a p y 
2 + 3 + 4 



1 



(5) 



B ± 



- B, f3 = A + 2B, y = -B and z = 7i(p LD ^f. 
This model and the linear one have been squared and fit to the 
data (Fig. 4). The results of this fitting up to the second lobe are 
summarized in Table 2. 




Projected baseline 



Projected baseline 



Fig. 4. Squared visibility data fitted by the best limb darkening model 
(dotted line). To improve the modeling, we should reproduce the excess 
of contrast at high frequencies. The squared visibility of a UD is drawn 
in full line for comparison. 



The third and fourth visibility lobes of Betelgeuse clearly 
contain high spatial frequency information more complex than 
just limb darkening and, up to this point, have been excluded 
from the fits. The linear limb darkening coefficient A=0.43 for 
a diameter of 44.28 ± 0.15 mas can be compared with the esti- 
mation by Perrin et al. (2004) in K band who found A=0.09 for 
a smaller diameter of -43.6 mas. Those measurements obtained 
in two different bands do not correspond to the same extension 
of the photosphere. That could explain the difference in diame- 
ter. On the theoretical side, Manduca (1979) measured A=0.52 
at 1.2 fim by fitting a simulated intensity profile of Betelgeuse 
with the following parameters: T e ff = 3750 K, log g = 1.5 and a 
solar composition. More recently van Hamme (1993) predicted 
A=0.49 for a limb darkening in J band with T e ff = 3500 K, log 
g = 0.5 and a solar composition. Given that limb darkening is 
probably smaller in the H band than in the J band, the agreement 
with our results is quite good. 

4.2. The MARCS modeling of the photosphere 

To model the photosphere of this supergiant, we use the 
SOSMARCS code, version May 1998, described in Verhoelst 
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et al. (2006) as developed by Gustafsson et al. (1975) at Uppsala 
University, Sweden (Plez et al., 1992, 1993). This code was 
specifically developed for the modeling of cool evolved stars, i.e. 
a lot of effort was put into the molecular opacities, and they allow 
for the computation of radiative transfer in a spherical geometry. 
A spectrum derived from this modeling is shown in Fig. 5. In 
that figure we note that the H narrow-band filter clearly does not 
see much water. Also CO and SiO molecules have been found 
in the MolSpheres around evolved stars, but neither molecule 
shows strong opacity in the H band. More interesting is the sug- 
gested presence of CN in a plume originating from the surface of 
Betelgeuse (Kervella et al., 2009). CN was not detected in these 
latter H-band data even if it does have an absorption band with 
some opacity across these wavelengths (see Fig. 9 in that paper). 
As their dynamic range was much higher than the one obtained 
in the observations analyzed here, CN is therefore not a concern 
for the results presented in this paper. From the brightness dis- 
tributions produced by this model, we can compute H-band vis- 
ibilities, to be compared with the IOTA observations. The result 
is shown in Fig. 6. 




1.5 2.0 2.5 

Wavelength [yu.m] 



Fig. 5. The photospheric spectrum (MARCS model) is plotted in black. 
The photosphere seen through the water layer suggested by Verhoelst 
et al. (2006) is in green. The optical depth of that layer (times 10 4 for 
visibility) is represented by the vertical lines, and the dashed line shows 
that the H-band filter is unaffected by circumstellar water opacity. 

The fit of the MARCS model leads to a photospheric limb- 
darkened diameter of 46.74 ± 0.12 mas (see Table 2). The re- 
solved flux ratio is estimated to ~5%. A Rosseland diame- 
ter (given by the distance from the star's center at which the 
Rosseland optical depth equals unity, tr oss = 1) of 45.03 ± 
0.12 mas can be deduced. 




10 20 




10 20 30 



Projected baseline 

Fig. 6. Upper panel: MARCS H filter squared visibility (in solid line) 
fitted to the data (dots) with diameter and resolved flux ratio as free vari- 
ables. Lower panel: residuals of the fit in the upper panel. We clearly 
see that this MARCS modeling is convincing up to the third lobe. 

ture at the surface of the star. It is a radiometric quantity that 
characterizes the mean surface brightness of the star. The spots 
imaged with IONIC on IOTA (Sect. 5.3) induce local variations 
of the surface brightness. It is interesting to compute how the 
effective temperature is sensitive to these spots. From our mea- 
surements, we deduce that the spots contribute to 2% of the ef- 
fective temperature, the same order as the uncertainty on the 
effective temperature. Current effective temperature measure- 
ments are therefore not accurate enough to detect spot fluctu- 
ations. 

4.4. Discussion 



Model 


0LD 


Limb darkening 


Reduced 






parameters ) 


x 2 


Linear 


44.28 + 0.15 


0.43 + 0.03 


11 


Quadratic 


44.31 + 0.12 


0.42 + 0.03 et 0.02 + 0.03 


11 


MARCS 


46.74 + 0.12 




12 



Table 2. Results of fitting of the limb darkening models. 



4.3. Effective temperature 

The effective temperature computed in Perrin et al. (2004) can 
be updated with these new diameter measurements. Following 
Dyck et al. (1992), a visual extinction of 0.5 was adopted in this 
calculation. We used the bolometric flux reported in Perrin et al. 
(2004): 111.67 ± 6.49 10" 13 Wcm -2 . Using these values, an aver- 
age diameter of 44.28 ± 0.15 mas for the limb darkening model 
and a Rosseland diameter of 45.03 ± 0.1 2 mas for the MARCS 
model lead to effective temperatures of3612±53K and 3587 ± 
53 K respectively. We propose the updated slightly lower value 
of 3600 ± 66 K for the effective temperature of Betelgeuse. The 
effective temperature has no meaning about the local tempera- 



The residuals show that there is a discrepancy at small and 
large baselines. At short baselines, we see that large residuals 
scatter in a symmetric manner around zero. This variation of 
visibility at these baselines is likely due to the effect of convec- 
tion on the diameter measurement at different position angles. 
The impact of convection on our interferometric measurements 
is being investigated by some co-authors of the present paper 
(Chiavassa et al., 2009b). It is remarkable that visibility data can 
be very well fit by a limb-darkened disk and that no MOLsphere 
is required, consistent with earlier studies in K (Perrin et al., 
2004). At large baselines, we see once again that the model can- 
not reproduce the observed high levels of visibility. 
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The use of a limb-darkening disk model greatly decreases the 
reduced^ 2 compared to a uniform disk model. A simple linear 
law seems to model this effect up to the third lobe. Beyond that 
limit, the data show a great complexity which is not reproduced 
by the simpler models. Contrary to AMBER K-band observa- 
tions made in January 2008 and recently published in Ohnaka 
et al. (2009) where the data points are scattered around a limb 
darkening model, the present observations reveal significant de- 
parture from symmetry. 

All limb darkening models that we used (linear, quadratic 
and MARCS model) hint at the presence of additional asym- 
metric structures suggested by the excess of contrast at high fre- 
quency and are beyond the scope of the limb-darkening models. 
In order to investigate this effect, we analyzed the closure phase 
measurements to determine the nature of asymmetric structures. 

5. Image reconstruction and asymmetries 

The dataset comprises 131 closure phase measurements. The 
first thing to note is that closure phases are not null over the 
whole dataset showing departure from point symmetry over 
Betelgeuse's stellar disk. Structures appear in closure phases 
when plotted versus triplet maximum projected baseline (Fig. 7). 



In the present work, we chose as an a priori object the best 
limb-darkened fit. Our prior is thus a purely symmetric limb 
darkened disk of 44.28 mas diameter surrounded by an environ- 
ment as found in section 4. In order to compare the two methods 
we made sure to use the same field of view, same sampling and 
same a priori image. Both use a quadratic regularization mean- 
ing that strong intensity gradients between the a priori image and 
the reconstructed image are quadratically discriminated against. 

5.1. MIRA 

MIRA (Multi-aperture Image Reconstruction Algorithm 
Thiebaut, 2008) is a method that iteratively minimizes a crite- 
rion to reconstruct the image. To minimize the criterion, MIRA 
uses the optimization method VMLMB (Thiebaut, 2002), a 
limited variable metric algorithm which accounts for parameter 
bounds. This last feature is used to enforce positivity of the 
solution. To avoid falling in local minima, the iteration starts 
with a strong regularization to the a priori image. Then the 
weight of the regularization decreases to favor the data fit (i.e. 
fi decreases). In terms of reduced the best reconstruction 
achieved gave a x 2 value of 5.7. The reconstructed image is 
shown in Fig. 8. 




-iq I- -I L_, i L i , i , i , i l: 

13.5 14.0 14.5 15.0 15.5 19.4 19.6 

Maximum projected baseline of a triplet (meters) Maximum projected baseline of a triplet (meters) 

Fig. 7. Closure phase measurements plotted versus two ranges of the 
triplet maximum projected baseline. 

With the lack of spatial frequency coverage of near-infrared 
interferometric observations, parametric modeling with simple 
geometrical descriptions turns out to be inefficient. Even with a 
few hundred visibility and closure phase measurements, it is dif- 
ficult to distinguish basic models because of the presence of a 
great number of local minima. Thus the parametric approach is 
not a good way to develop a complex model. In the last decade, 
various teams have therefore undertaken the creation of image 
reconstruction algorithms adapted to interferometric data in the 
near infrared. In this paper we used two of them: MIRA and 
WISARD. Both are presented and compared in Le Besnerais 
(2008) and in Lacour et al. (2008). The major advantage of such 
methods is that they can reveal structures by reconstructing an 
image which fits the data and which is also close to an a priori 
image of the object. The balance between the regularization of 
the reconstructed image to the a priori image and the data fitting 
is allowed by setting a hyperparameteryu: 

J total ~ J data ~^ I^J apriori (6) 

where J data* Japrwri are the criteria that quantify how the recon- 
structed image fits the data and how it is close to the a priori 
image. J ma \ is the final quantity to minimize. 



5.2. WISARD 

WISARD (Weak-phase Interferometric Sample Alternating 
Reconstruction Device, Meimon , 2005; Meimon et al., 2009) 
has a more explicit approach to the missing phase information is- 
sue. Using a self calibration method inspired from radio interfer- 
ometry (Corn well & Wilkinson, 1981), it aims at reconstructing 
an image and also finding the missing phase data. WISARD uses 
an appropriate approximation of the specific noise encountered 
in infrared interferometric data, which favors the convergence 
toward a global minimum - see Meimon et al. (2005) for details. 

For the best reconstruction obtained with WISARD, the min- 
imum value of this criterion corresponded to a reduced^ 2 of 4.7. 
The reconstructed image is presented in Fig. 8. 

5.3. Analysis of the reconstructed images 

The first thing to note is the impressive similarity between the 
two images given the complexity of the data. The two algorithms 
both reveal two spots (named Tl and T2) of unequal brightness 
and which are roughly located at the same positions (Table 3). 
Apart from Tl and T2, the reconstructed images look very sim- 
ilar. The fact that two independent methods give similar results 
expresses their consistency. 

We tried different regularization types and saw that the re- 
sulting images were not different from the quadratically regular- 
ized ones. We believe the reality and the location of asymmetries 
on the stellar disc of Betelgeuse can be trusted. The interferom- 
eter image of a point source was also reconstructed (Fig. 8): the 
asymmetries imaged at the surface of Betelgeuse are not arte- 
facts caused by the sparse uv plane coverage. A visibility plot 
comparing the data and the reconstructed visibilities of MIRA is 
presented in Fig. 9. The fit of the V 2 data is clearly very good 
and the best of those obtained in this paper. 

In Table 3, we see that Tl is bigger than T2 and brighter 
by a factor of 2. Tl seems more elongated in the northwest- 
ern direction than T2 though the inclination changes between 
the two images. The limiting resolution we achieved during our 
observations is about 1 1 mas. Separated by a distance of about 
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Parameter 


1V11KA- 1 1 


Wlo/VKJJ- 1 1 


1V11KA- 1 Z 


Wlo/VKJJ- 1 Z 


Iv Iv /A - L U 1 1 


z,.j in I.J 




-L.D It l.J 


f\ + 1 s 

-u.u It L.J 


RDEC-cen 


7.5 + 1.5 


7.1 + 1.5 


-4.5 + 1.5 


-4.7 + 1.5 


FWHM-RA 


10.32 + 2 


11.08 + 2 


6.86 + 2 


8.89 + 2 


FWHM-DEC 


5.80 + 2 


6.38 + 2 


4.32 + 2 


5.65 + 2 


Amplitude 


8 + 1 


9 + 1 


4+ 1 


5 + 1 


(% the max) 











Table 3. Parameters of a Gaussian fit to the two bright spots revealed by MIRA and WISARD. Tl corresponds to the northern spot, T2 to the 
southern one. 



12.5 mas, the asymmetry created by an unequally bright pair 
of spots such as Tl and T2 was resolved by the interferometer. 
With an angular size of about 1 1 mas, Tl is resolved and we can 
trust the presence of this spot at this location on the stellar sur- 
face. Tl might also be a combination of smaller scale structures. 
However it is not the case for T2 which is smaller than the max- 
imum resolution we got. We therefore cannot conclude on the 
size that T2 may intrinsically have. Betelgeuse is certainly more 
complex than it appears on these images. Particularly the pres- 
ence of structure whose size is smaller than 10 mas is very likely 
when looking at these images. However, the respective reduced 
X 2 values of MIRA and WISARD are low enough that it can 
be concluded that the reconstructed images (respectively 5.7 and 
4.7) represent a fairly convincing modeling of the star surface 
given the accuracy, the complexity and the amount of data. 

5.4. Parametric fit of two spots 

Image reconstruction has unveiled unforeseen structures with the 
use of an a priori image and some regularization conditions. 
Nevertheless, fitting a parametric model is an important task 
since it models the object with a limited number of parameters 
without the influence of the regularization. As the two different 
regularization methods point towards the existence of two bright 
structures near the center of the star, we used this information 
to simultaneously model visibilities and closure phases. We thus 
modeled two spots as uniform disks plus the limb-darkened disk 
and environment derived in section 4. The diameter of the north- 
ern spot was set to the value found in the previous section. The 
free parameters were the positions and relative fluxes of the two 
spots and the diameter of the southern spot. Results are shown 
in Fig. 10 and in Table 3. 

The best fit gives a reduced^ 2 of 15 for the spot positions 
given in Table 3. The flux ratio and diameters corresponding to 
the best fit are listed in Table 4. 



Parameter 


Param-Tl 


Param-T2 


RRA-cen 


6.5 + 1.5 


-6.2 + 1.5 


RDEC-cen 


6.8 + 1.5 


-8.0+ 1.5 


UD diameter in mas 


11 


0.02 + 1 


Amplitude 


2+ 1 


2 + 1 


(% the max) 







Table 4. Parameters of the 2 spots found by the parametric fit. Spots are 
modeled as uniform disks. 



The third visibility lobe fit is clearly improved with a two- 
spot model compared to circularly symmetric models. It is con- 
firmed that the southern spot is not resolved although the flux 
ratios found by this method do not match those derived from the 
regularized images. Contrary to the case of the regularized im- 
ages, the fourth lobe is not well reproduced by the parametric fit, 



suggesting that even smaller scale bright structures are needed to 
interpret these data. This may partially explain the discrepancies 
with the regularized images which better fit the data. More gen- 
erally, all the structures disclosed by reconstruction algorithms 
are very difficult to model and fit in a parametric way because 
of the number of free parameters required. Although the 2-spot 
model catches most of Betelgeuse surface brightness complex- 
ity and can be considered a fair description of the object spatial 
brightness distribution, going further requires even smaller and 
fainter details (as images in Fig. 8 show) thus probably showing 
the complexity of the star surface. 

6. Conclusion 

In this paper, we report on high accuracy interferometric mea- 
surements of Betelgeuse obtained with IOTA/IONIC up to a res- 
olution of ~1 1 mas. The precision on squared visibilities enabled 
by the use of guided optics was better than a percent. The com- 
plexity revealed by the data led us to use different methods and 
models which gave important results. We: 

- determine the limb-darkened diameter of Betelgeuse using 
parametric and MARCS models as well as a new effective 
temperature, 

- detect a fully resolved circumstellar environment as bright 
as 4% of the total flux. It is very likely that the silicate dust 
shell is a significant contributor to this environment, and 

- image two spots at the surface of the star. One of the spots 
is unresolved and the other one is as large as one half of the 
stellar radius. 

No MOLsphere has been detected in agreement with previ- 
ous work at near-infrared wavelengths. We thus find a compat- 
ibility with the predictions published in Schwarzschild (1975) 
about the existence of large structures on the surface of 
Betelgeuse due to convection. As Kiss et al. (2009) observed for 
other red supergiants, we note that these structures are responsi- 
ble for up to -10 percent of the total flux. 

If one assumes that Betelgeuse and the spot Tl radiate like 
blackbodies, an effective temperature of 3600 K for Betelgeuse 
leads to a spot Tl effective temperature of 4125 K (11 mas 
diameter and a flux ratio of 0.085). This is compatible with 
Schwarzschild (1975) who predicted up to 1000 K temperature 
difference between the hottest and coolest parts of a convective 
surface. If these spotty structures are similar to those observed 
on the Sun, then they could be faculae which are bright zones 
surrounding dark convection cells. Their presence on Betelgeuse 
would show that a magnetic field plays an important role as it is 
the case for the Sun. 

To go further, it seems that the analysis of this dataset could 
be improved by injecting small scale (< 10 mas) structures 
in the modeling. This can be done by comparing our data to 
simulated visibilities and closure phases generated from RSG 
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Fig. 8. Upper window: Contour image reconstruction from MIRA. 
Middle window: Contour image reconstruction from WISARD of 
Betelgeuse in a 60 mas field. Both images were reconstructed with the 
same a priori object and the same type of regularization. Lower win- 
dow: Image of the dirty beam, i.e., the interferometer's response to a 
point source, given our uv plane coverage. 



hydro-radiative simulations (such as presented in Chiavassa et 
aL, 2009a). This would undoubtedly bring new clues for the con- 
vective nature of the structures found in the present work. The 
convective motions could play an important role in the mass-loss 
process operating in RSGs, because of the strong decrease in ef- 
fective gravity due to turbulent pressure (Josselin & Plez, 2007). 
The bright plume found by Kervella et al. (2009) might be the 
first direct visible evidence of a mass loss mechanism in RSGs. 
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Fig. 9. Reconstructed squared visibilities from MIRA compared to the 
observed data. Only the second, the third and the fourth lobe are plotted 
here. We note that the third and fourth lobes are now well reproduced 
thanks to the contribution of two bright spots. The WISARD recon- 
structed visibilities are very similar to the MIRA ones and are not plot- 
ted here for conciseness. 



Its formation could be explained by large upwards moving con- 
vective cell or by rotation, through the presence of a hot spot at 
the location of the polar cap of the star. The possibility that this 
plume is connected to a spot on Betelgeuse should be investi- 
gated with interferometric techniques in the near future. These 
perspectives hold the promise of determining a first estimation 
of the convection scale in RSGs which would represent a new 
step in the understanding of the mass loss enigma in RSGs. 
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